The excitation of rogue waves in a variable medium: an experimental study on the 

interaction of water waves and currents 
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We show experimentally that a stable wave propagating into a region characterized by an opposite 
current may become modulationaly unstable. Experiments have been performed in two indepen- 
dent wave tank facilities; both of them are equipped with a wavemaker and a pump for generating 
a current propagating in the opposite direction with respect to the waves. The experimental results 
support a recent conjecture based on a current-modified Nonlinear Schrodinger equation which es- 
tablishes that rogue waves can be triggered by non-homogeneous current characterized by a negative 
horizontal velocity gradient. 

PACS numbers: 



Ocean waves are characterized by a statistically small 
steepness and often (but not always, see for example [H) a 
weakly nonlinear approach is sufficient to capture some of 
the intriguing aspects hidden in the fully nonlinear prim- 
itive equations. This weakly nonlinear approach is also 
shared by other fields of physics such as for example non- 
linear optics Q and plasma physics where small pa- 
rameters can be individuated and asymptotic expansions 
can be used to simplify the original equations. If the con- 
sidered physical process is not only weakly nonlinear but 
also narrow-banded then the lion's share is played by the 
Nonlinear Schrodinger equation (NLS). Being an exactly 
intcgrable equation via the Inverse Scattering Transform 
[H, bizarre analytical solutions have been found in the 
past: besides traveling waves, breathers or multi-breather 
solutions have been found [5|-|7| and observed in hydro- 
dynamic [1, Q , nonlinear optics fiol . [Tl| and plasma pj| 
experiments. Starting from fl3l Il4j . such solutions have 
been considered as prototypes of rogue waves. The earl 
stages of the so called Akhmediev breather solution 
describes the exponential growth of a slightly perturbed 
plane waves, i.e., it corresponds to the classical modula- 
tional instability process [15( . For water waves in infinite 
water depth, the instability is active when e N > 1/ v2 5 
where e = fcrjflo is the initial steepness of the plane wave, 
kg the wavenumber, ao its amplitude and TV = coo/Afl is 
the number of waves under the modulation; loq is the an- 
gular frequency corresponding to the carrier wavenumber 
/co and AQ the angular frequency of the modulation. 

The whole picture is by now pretty well understood 
and relies on the fact that the medium in which waves 
propagate is homogeneous. In terms of the NLS equa- 
tion this means that the coefficients of the dispersive 
and nonlinear terms do not depend on the spatial co- 



ordinates. Much more complicated and intriguing is the 
case in which the medium changes its properties along 
the direction of propagation of the waves. This situa- 
tion is much more difficult to treat in terms of simplified 
models because it turns out that in general the resulting 
modified NLS does not share the property of intcgrability 
as the standard NLS and analytical breathers solutions 
can be found only in special cases (see some examples in 

EMU). 

In the oceanographic context, the non- homogeneity of 
the medium is mostly due to currents or bottom topogra- 
phy. In this Letter we will consider, from an experimental 
point of view, the interaction of waves and currents and 
the consequent formation of rogue waves. Wave-current 
interaction and its effect on wave instability have been 
the subject of a number studies over the past decades 
[l9l - l26l ]. Only recently, however, has the rogue wave for- 
mation process induced by an opposite current gradient 
been studied numerically on the basis of the modified 
NLS equation [27M30j . Specifically, [28[ applied a trans- 
formation to the one dimensional current-modified NLS 
equation derived in (27} to obtain the following NLS-type 
of equation: 



dB k d 2 B 



dx 
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+ ik^exp(-AU/c g )\B\ 2 B = 0, (1) 



where c g is the group velocity, AU = U(x) — U(0), with 
U (x) the velocity of the current at position x and U(0) 
is the current at x = 0. For simplicity, we will consider 
the physical case of a wave generated in a region of zero 
current, U(0) = 0, that enters into a region where an 
opposite current starts increasing its speed (in absolute 
value) and then adjusts to some constant value Uq- Note 
that in this case the coefficient of the nonlinear term of 
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equation (fT]) increases as the waves enter into the current 
up to a certain value and then remains constant. The net 
effect is therefore an increase of the nonlinearity of the 
system. Numerical simulations of the current modified 
NLS equation presented in [28| showed that an envelope 
of an initially stable wave train becomes unstable after 
entering in the current region. As a result, the maximum 
amplitude shows a growing trend for increasing the ratio 
Uo/cg, corroborating the idea that an originally stable 
plane wave is transformed into a breather by the presence 
of a current . 

In [3(| it was noted that the modified NLS equation 
proposed in [27l and used in [28| does not preserve wave 
action (see |2l[ ). To the lowest order, conservation of 
wave action can be accounted for in equation ([T]) by sim- 
ply multiplying the ratio Uo/c g by a factor 2. Based on 
this modification, a prediction for the maximum wave 
amplitude obeys the following equation: 



E 



= 1 




CXP (Up/Cg) 

V2eN 



(2) 



where A max is the maximum wave amplitude achieved in 
the region of constant current and \f~E is standard devia- 
tion of the wave envelope once the current has reached its 
maximum constant value. In (30j a derivation of a modi- 
fied NLS equation based on an Hamiltonian formulation 
of surface gravity waves has been performed. A similar 
prediction as the one in ^ has been proposed and takes 
the following form: 
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Note that, to the lowest order, an expression equivalent 
to ([3]) can be derived from equation ([1]) by multiplying 
the ratio Uo/c g by a factor 3. 

In this Letter we present two independent sets of lab- 
oratory experiments that were conducted in the wave 
flume of Plymouth University and in the narrow direc- 
tional wave basin at the Ocean Engineering Tank of the 
University of Tokyo. The wave flume at Plymouth Uni- 
versity is 35 m long and 0.6 to wide with a uniform water 
depth of 0.75 m. The facility is equipped with a piston 
wavemaker with active force absorption at one side and 
a passive absorber panel at the other end; only unidirec- 
tional propagation is allowed. The flume is also equipped 
with a pump for the generation of a background current 
up to 0.5 m/s, which can follow or oppose the wave direc- 
tion of propagation (only an opposing current was used 
for the present study though). One of the inlet/outlet 
is located nearby the absorber, while the other is at a 
distance of about 2.5 to from the wavemaker. This par- 
ticular configuration allows waves to be generated outside 
the current field and propagate for a few wavelengths be- 
fore encountering a current gradient. The wave field was 
monitored with 10 wave probes equally spaced along the 
tank, while the velocity field was monitored with two 



Acoustic Doppler Velocimctcrs (ADV) properly seeded. 
A survey of the current revealed a fairly uniform flow 
both in space and time. An example of longitudinal and 
vertical profiles of the average horizontal velocity is pre- 
sented in Figs. QJi aiidQJ). 

The Ocean Engineering Tank of the Institute of Indus- 
trial Science, University of Tokyo (Kinoshita Laboratory 
and Rheem Laboratory), is 10 m wide, 50 to long and 5 m 
deep. It is equipped with a multidirectional wavemaker 
with 32 triangular plungers (0.31 m wide) [3l|. A sloping 
beach is deployed opposite the wavemaker to absorb the 
wave energy. The tank is also equipped with a pump (lo- 
cated beside the basin) for the generation of background 
currents up to 0.4 m/s; the stream can follow or oppose 
the waves. One of the inlet/outlet is located below the 
beach, while a second is located below the wavemaker 
(approximately 2 m below the water level). Note that no 
modification of the cross section was performed to locally 
modify the velocity field. Wave probes were deployed 
along the tank at a distance of 2.5 to from the sidcwall 
and arranged at 5 to intervals to monitor the evolution of 
wave trains (a six-probe array was also deployed to moni- 
tor directional properties) ; for consistency with the wave 
flume experiments, only probes within 25 to from the 
wavemaker were considered, though. Two electromag- 
netic velocimeters were used to survey the current. In- 
struments were deployed at several locations in the tank 
at a depth of 0.2 to; a vertical profile was measured at 
about 10 to from the wavemaker. Instantaneous mea- 
surements of horizontal velocity revealed that the flow 
had substantial spatial and temporal variations, with a 
dominant oscillation period of approximately 150 s. Av- 
erage values are presented in Figs. [TJ; and [Hi. For these 
tests, the standard deviation was about 25% of the mean 
over the entire time series. As the flow's outlet is located 
just below the wave generator, the velocity is approxi- 
mately zero at a distance of about 0.2 to from the wave- 
maker, while the flow is at regime at a distance of 5 to 
from the wavemaker. Waves are therefore generated in 
a condition of (almost) still water and undergo a current 
gradient about 1 to after being generated. Farther from 
the wavemaker, between 5 and 30 to from the genera- 
tor, the current still shows a weak gradient, which may 
slightly affect the wave field. Although the average hor- 
izontal velocity weakly decreased with the water depth, 
the vertical profile remained fairly uniform over a depth 
of about 1.5 to. The survey of the current field also indi- 
cated that the current runs faster on one side (along the 
wave probes), while it is slower on the other (see Fig. Q}:, 
for example). It is important to mention that the cross- 
tank gradient refracts the wave field towards the sidcwall 
(probes side). As such, this may cause linear directional 
focusing, steepening the wave profile and hence affecting 
wave dynamics. 

A number of tests characterized by different values of 
the modulation frequency and current velocity have been 
carried out in both facilities. For all tests, the initial 
signal at the wavemaker consisted of a three-component 



3 



E 

~-0.05 
-0.1 
-0.15 
-0.2 



B 









Q. 


-0.2 


CD 


Q 






-0.4 




-0.6 



(a) 



E 

3-0.05 
-0.1^ 



-0.15^ 
-0.2 



10 20 30 

Distance from the wavemaker (m) 

Oft # at 2.5 m from side wall i 
□ at 5.0 m from side wall " 
Aat 7.5 m from side wall 

6a - A.. A A A 

w » y w :S 



0i 



A 



(c) 

10 20 ^30 

Distance from the wavemaker (m) 



-1.5 



(b) 

-0.2 -0.1 

U (m/s) 

• 

(d) 



40 



-0.2 -0.1 

U (m/s) 



FIG. 1: Average horizontal velocity for a current field op- 
posing wave propagation: longitudinal (a) and vertical (b) 
profiles in the wave flume of Plymouth University; and lon- 
gitudinal (c) and vertical (d) profiles in the directional wave 
basin of the University of Tokyo. The longitudinal profiles 
were measured at a depth of 0.3 m at the University of Tokyo 
and 0.08 m at Plymouth University. 



system: a carrier wave of period Tq = 0.8 s (wavelength 
Ao = 27r//co — 1 m) and two side bands with amplitudes 
b± equal to 0.25 times the amplitude a c of the carrier 
waves. Considering the water depth of 0.75 m in the wave 
flume and 5 m in the directional basin, experiments were 
performed under deep water conditions (fcn/i > 4). As 
the effect of the current is to steepen the wave profile, a 
small initial steepness was selected in order to avoid wave 
breaking. The tests presented here were carried out by 
selecting the wave amplitude of the carrier wave a c in 
such a way that the wave steepness was fcnQo = 0.063 
with a§ = a 2 + b\ + b 2 _. The frequency of the distur- 
bances was chosen to force the number of waves under 
the perturbation N = ujq/AQ (with ujq being the angular 
frequency of the carrier waves) to be equal to 11. Un- 
der these circumstances, the perturbation frequency lays 
just outside the NLS-based instability region, i.e. waves 
are stable (eN = k Q a N = 0.69 < 1/V2). These pack- 
ets were tested against different opposing current speeds. 
Velocities ranged from m/s to -0.30 m/s with step of 
-0.02 m/s. All tests were run for a time period of 10 
minutes. Considering the variability of the current in the 
directional basin, this ensured enough data to perform a 
standard statistical analysis of the observations. 

Typical time series of the water surface elevation are 



presented in Fig. [2] for U~o/c g = and Uo/c g = —0.1. 
What is clear from the figure is that, while the modula- 
tion does not grow in absence of current (top panels) , the 
wave packets undergo a modulationally unstable process 
(nonlinear focussing) and wave amplification in the pres- 
ence of a current. Interestingly enough, amplification is 
less sharp in the wave flume due to a more regular current 
field and the absence of contaminating three-dimensional 
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FIG. 2: Evolution of surface elevation: samples from the wave 
flume of Plymouth University (right panels) and the direc- 
tional wave basin of the University of Tokyo (left panels). 



effects. 

A standard zero-crossing procedure was applied to ex- 
tract the maximum amplitude at each probe. Because of 
the temporal variability of the current in the directional 
basin, the analysis was performed on segments of three 
consecutive wave groups, i.e. time window r = 26.4s, 
where the current was nearly steady (velocity fluctua- 
tions were confined within a range of ±0.01 m/s, namely 
one order of magnitude smaller than the average). As 
the prediction in equation <j2j) only includes the con- 
tribution of free wave modes, frequencies greater than 
1.5 w c and smaller than 0.5w c were removed to filter 
out bound modes. The amplitude was then normalised 
by E 1 ' 2 = [(1/t) J t \A\ 2 dt}^ 2 , where A is the wave 
envelope of the concurrent segment, to eliminate the 
current-induced increase of wave amplitude. An average 
normalised maximum amplitude and standard deviation 
were calculated over the entire time series. The maxi- 
mum wave amplitude is presented as a function of Uq /c g 
in Figure [3] together with equation @ (solid line) and 
the prediction model in ^ (dashed line); error bands 
equivalent to the 68% confidence interval (one standard 
deviation) are also shown; owing to the stable current 
field, the error band for flume experiments is notably 
smaller than the one detected in the directional basin. 
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FIG. 3: Normalised maximum amplitude as a function of 
Uo/c g : equation |[2j, solid line; equation ([2]) modified by in- 
cluding the conservation of wave action, dashed line. 

Qualitatively, both tests are in good agreement with the- 
ory, substantiating that an adverse current gradient trig- 
gers moduational instability processes. Quantitatively, 
however, observations are notably overestimated by the 
model in ([3]), especially for strong current fields. Equa- 
tion on the contrary, produces a satisfactory approx- 
imation of the records, particularly for the flume exper- 
iments. It is important to remark, in this regard, that 
three dimensional effects induced by refraction enhance 
the breaking probability in the directional wave basin 
and hence limit the maximum wave amplitude also for 
relatively mild currents. Amplitude growth ceases, in 
fact, for \Uo/c g \ > 0.2 as individual waves systemati- 
cally reach their limiting steepness [32^ and break. In 
the wave flume, three-dimensional effects are suppressed, 
making breaking dissipation less likely and thus allowing 
waves to develop the maximum amplification. This justi- 
fies the good agreement with theory until |[/o/c 9 | ~ 0.4; 
beyond this threshold waves reach the limiting steepness 
and break. 

In conclusion, the conjecture on the generation of rogue 
waves in opposite current using a current-modified cu- 
bic NLS equation has been here confirmed by two inde- 
pendent sets of laboratory experiments, which were car- 
ried out in a wave flume and a narrow directional wave 
basin. Observations corroborate that the excitation of 
the modulation and the concurrent intensification of the 
maximum wave growth is a function of the adverse cur- 
rent Ug/cg and is consistent with theoretical predictions 
in equation @, but overestimated by the model in ([3]). 
Despite possible shortcomings of the theory and uncer- 



tainties in the experimental conditions, especially related 
to the irregularities of the current field, the present study 
clearly shows evidence that an opposing current field can 
destabilize an otherwise stable wave packet. The result 
is the formation of a rogue wave, whose amplitude de- 
pends on the ratio of current velocity to group velocity. 
Note that the essence of the theoretical predictions in [28[ 
and [30l ] is to account for the impact of the current shear 
only to modify the basic parameters that control the dy- 
namics of the wave train. Once the basic parameters 
have changed, the instability process follows its nominal 
evolution pattern (namely, the nonlinear stages of modu- 
lational instability). It is important to mention that this 
is the case for our experiments because the spatial scale 
of the current field variation was smaller than the spatial 
scale of the nonlinear evolution. This is not necessar- 
ily the case in reality and further research is needed for 
those cases when the two scales are of the same order. 
From a physical point of view, the process that we have 
described may take place in nature when a modulation- 
ally stable swell (waves propagating without the forcing 
of the wind), which is characterized by a narrow spec- 
trum (both in direction and frequency) enters a region 
of an opposite variable current. Besides the linear effect 
of refraction which could generate linear focussing, the 
current gradient can destabilise the wave packets leading 
to a nonlinear focussing effect and the formation of rogue 
waves. Indeed, as mentioned in [28| . such currents may 
reach velocities up to 1.5 m/s, and for a group velocity 
corresponding to waves of period equal to 10 s (a typ- 
ical condition during storms), the ratio Uo/c g is of the 
order of 0.2, large enough to trigger a dangerous rogue 
wave. We believe that the mechanism we have observed 
is universal and can be reproduced in all fields of physics 
where an NLS equation with variable coefficient of the 
nonlinear term can be written. 
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